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The asymmetric simple exclusion process with open boundaries, which is a very simple model of 
out-of-equilibrium statistical physics, is known to be integrable. In particular, its spectrum can be 
described in terms of Bethe roots. The large deviation function of the current can be obtained as 
well by diagonalizing a modified transition matrix, that is still integrable: the spectrum of this new 
matrix can be also described in terms of Bethe roots for special values of the parameters. However, 
' due to the algebraic framework used to write the Bethe equations in the previous works, the nature 

of the excitations and the full structure of the eigenvectors were still unknown. This paper explains 
C , why the eigenvectors of the modified transition matrix are physically relevant, gives an explicit 

^ ' expression for the eigenvectors and applies it to the study of atypical currents. It also shows how 

the coordinate Bethe Ansatz developped for the excitations leads to a simple derivation of the Bethe 
CO , equations and of the validity conditions of this Ansatz. All the results obtained by de Gier and 

CN ■ Essler are recovered and the approach gives a physical interpretation of the exceptional points The 

overlap of this approach with other tools such as the matrix Ansatz is also discussed. The method 
that is presented here may be not specific to the asymmetric exclusion process and may be applied 
^ . to other models with open boundaries to find similar exceptional points. 
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■ I. INTRODUCTION 

C/2 

Out-of-equilibrium statistical physics has been an active field of research for more than twenty years. Breaking 
the detailed balance property or imposing currents of particles or energy through a system prevents from applying 
I most of the standard techniques of equilibrium statistical mechanics. The determination of the complete distribution 
. of the currents, beyond the average values, is in general a hard problem since it depends on the whole history 

■ followed by the process. Many general symmetry properties [11, 0, 0] have been established for a wide class of models 
O and the macroscopic fluctuation theory 0, makes the computations of some macroscopic quantities possible for 

, ^/ diffusive systems. However, it is hard to compare these results with the ones that can be obtained from a microscopic 
\ approach since there are only a few examples of exactly solvable models. This comparison becomes necessary for 

■ far-from-equilibrium models for which no general macroscopic theory exists yet. 

^ ' Finding out-of-equilibrium models that can be exactly solved is often less easy than for equilibrium models. One 
00 of the reasons, which comes from integrability, is that the out-of-equilibrium character of the system comes from the 
boundary conditions imposed to the system. Putting the system in contact with reservoirs at different temperatures 
or different chemical potentials requires to have a description of the interaction between the system and each reservoir. 
From the point of view of integrability, the diagonalization of the Markov transition matrix of the system requires 
\ that both the dynamics of the bulk and the one of the contact with the reservoirs have to be integrable, while at 
' equilibrium only the bulk dynamics may be relevant. 

' The asymmetric simple exclusion process (ASEP) is one of the simplest out-of-equilibrium models that is integrable. 
• • The ASEP describes hard-core particles that diffuse on a one-dimensional lattice of L sites. There is at most one 
particle per site. Each particle tries to jump to the right (resp. left) neighboring site with rate p (resp. q): if the 
target site is empty then the jump is performed, if the site is occupied then the particle remains on the initial site 
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Figure 1: Transition rates in the asymmetric exclusion process for a system of size L. 
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(exclusion effect, see fig. [T]). The first and the last sites are in contact with reservoirs of particles that can add or 
remove particles in the system so that they tend to impose their own densities in the system. If the two densities are 
different and/or if the bias p — q is non-zero, then a current of particles flows through the system from one reservoir 
to the other. 

Many properties of this one-dimensional model of particle transport are now known. A first breakthrough was the 
introduction of the matrix ansatz Q to study the stationary measure of the ASEP with open boundaries. Although 
this technique has been fruitfully generalized to many other systems f?^ , it remains confined to the study of stationary 
measures and does not give any information about the relaxation times of the system nor the distribution of the 
current that flows through the system. The st udy of finite dimensional representations of the matrix algebra [1, Q 
and their interpretation in terms of shocks [Tol . [Til . [T^ give a new physical point of view on the matrix ansatz and is 
related to the approach I describe below. 

Besides the emergence of numerical methods to study the larg e deviation function of the current in the system 
[13I [13, [iB| , a second breakthrough has been obtained in [H, [ij HI] . The Markov transition matrix of the ASEP 
with reservoirs can be related to the Hamiltonian of the XXZ spin chain with non-hermitian coniplex boundary flelds 
and to the Temperley-Lieb algebra and many results obtained for the XXZ spin chain Elj US, US, UHl o^' 

the Temperley-Lieb algebra [2a . [27j can be extended to the ASEP. A major inconvenient of these methods is that 
they are often valid only at exceptional points in the parameter space and the validity condition did not have yet 
any simple physical interpretation in the case of the ASEP. Moreover these methods do not give the structure of the 
eigenvectors of the Markov transition matrix and thus prevent the computation of correlation functions. This paper 
presents a simple construction of the eigenvectors based on simple remarks on the expected properties of vacuum 
states and excitations. These eigenvectors take the form of a coordinate Bethe Ansatz. 

This paper is organized as follows: I introduce in sectionlTTl general properties of the ASEP and show how the current 
distribution can be obtained and how the coordinate Bethe Ansatz works for periodic boundary conditions. In section 
IlICi I explain the physical interpretation of the eigenvectors build in the next sections. Section |lll] is devoted to the 
study of the bulk dynamics and I show how to construct relevant vacuum states and excitations in presence of two 
different reservoirs at each end of the lattice. The structure of the excitations in part of the spectrum is similar to the 
shock structure described in [13) [O, [l^l > whom I refer for a physical discussion of the shocks, but I will not focus on 
it since the construction may be valid for many other systems with different interpretations of the excitations. Section 
IIVI focuses on the dynamics of an excitation that arrives on a boundary : the Bethe equations are derived and 
the interpretation of the boundary integrability is discussed. It also contains the complete formula for the proposed 
coordinate Bethe Ansatz. The last section IVl is an application of section ITVl and shows how the construction of these 
eigenvectors can be used to study and interpret the dynamics of the system conditioned to produce an atypical current 
as introduced in section HTCl 



II. GENERAL PROPERTIES 



A. Definition of the process and generating function of the current 



The asymmetric simple exclusion process corresponds to particles that diffuse on a lattice. We consider here the 
case of a one-dimensional lattice of length L. There is at most one particle per site. If the target site is empty, each 
particle can jump to the neighbouring site on the right with probability pdt and to the one on the left with probability 
qdt during an infinitesimal time dt. If the target site is occupied, the particle does not attempt to jump. An isolated 
particle performs thus a random walk. Through the bond between the site i and i -|- 1, the transition matrix in the 
basis of configurations (0^0^+1, Oili+i, liOi+i, lili+i) is thus given by : 

^0 ON 

^M+i - I n 7 i". n 1 ■ (II-l) 



q -p 
kO Oj 



One can notice that the number of particles is conserved by this bulk dynamics. In case of periodic boundary 
conditions, the site L is connected to the site 1 by a similar local matrix wl.i- The diagonalization of the Markov 
transition matrix using the coordinate Bethe ansatz is presented below in section IIIDI 

For open boundaries, each reservoir tends to impose its own density pa or ph on the neighboring site. The interactions 
with the reservoirs are described by boundary operators Bi and Bl acting on the sites 1 and L. If site 1 is occupied, 
the particle is removed with a rate 7; if it is empty, a particle is injected with rate a. In the same way, a particle is 
injected by the second reservoir on site L with rate S and is removed with rate /3 (see fig. [J). The two densities the 
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reservoirs try to impose are thus given by pa = + 7) and pi, = 5/{(3 + 5). In the bases (Oi, li) and (Ol, Iz,), the 
operators Bi and B]^ are described by the two matrices: 



B^ = \ : ' ■ (n-2a) 



—a 


7 


a 


-7 


-5 





= ( / _^ 1 ■ (n.2b) 



The Markov transition matrix for the open chain is thus given by: 

L-l 

W = Bi®7®^-i + ^/«*-i«)Wi,i+i®7®^-'-i+/®^-i®Bi, (II.3) 

i=l 

where / is the 2x2 identity matrix. If one introduces the vector \Pt) = X^c -Pt(C|Co)|C) where Pt{C\Co) is the 
probability that the system is in the configuration C at time t knowing that the system starts in configuration Co and 
the vectors \C) form a basis of a 2^-dimensional vector space, then one has the evolution 

^^\P,)=W\Pt). (II.4) 

The diagonalization of W gives the stationary measure and the relaxation times of the ASEP with open boundaries 
but not the full distribution of the current. Since the system is one-dimensional and the number of particles is 
conserved in the bulk, the current Q can be measured through any bond and will be defined here as the difference 
between the integer numbers of particles added on the first site by the left reservoir and of particles removed by the 
same reservoir. To obtain it, one has to consider the joint probability Pt{C,Q\Co) that the system is in configuration 
C at time t and that a total current Q has been counted between the initial time and t knowing that the system starts 
in configuration Cq. This probability evolves as: 

jPt{C,Q\C^) = J2 Wcc'Pt{C',Q- qcc'\Co) - I J2 ^c'c Pt{C,Q\Co), (II.5) 

C'jtC \C'#C / 

where qcc = ±1 or is the number of particle exchanged with the left reservoir during the change of configuration 
C C. The generating function of Q is defined as Pt{C,s\Co) = e*'^Pt(C, Q|Co) and satisfies the differential 
equation: 



d_ 
di 



Pt{C,s\Co) = J2 (Wcc'e''"''')Pt{C',s\Co) - [ ^ Wcc ] Pt{C,s\Co). (II.6) 

C'jiC \C'jtC J 



In this linear algebra formalism, the long time behavior of Pt{C, s\Co) is obtained by diagonalizing a modified matrix 
W obtained by replacing the boundary operator Bi by: 

This matrix W is not stochastic anymore and the long-time behavior of Pt{C, s\Co) is dominated by the first eigenvalue 
/Ui (s) of the matrix W, 

Pi(C,s|Co)(xe''iW*, (II.8) 

which corresponds to the large deviation behavior: 

Pt(C,g|Co)oce*/W/*), (II.9) 

where / and /ii are related by the Legendre transformation /ii(s) = maxj{f{j)+sj). Most of the interesting properties 
are thus contained in the modified matrix W and a way of diagonalizing it with Bethe ansatz methods is presented 
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in the next sections. The integrability of W is related to the Temperley-Lieb algebra formed by the operators Wi^i+i, 
Bi and Bl' 

=-(p + <7K.+i: = -(a + 7)^1, BI = -{P + 5)Bl, (Il.lOa) 

Wi,i+lWi+lA+2W^^^+l = (pq)Wi^i+i, Wi+l,j+2Wi,i+i Wi+i,,+2 = (TO)«^i+l,J+2, (II. 10b) 

Wi,2-Biwi,2 = -(ag +P7)wi,2, wl-ixBlwl-i,l ^ ~{Pq + pS)wl-i,l- (II.lOc) 

Many interesting properties of this algebra have been used in the context of integrable systems [H, and loop 
models [l^ . For example, the Yang-Baxter equations can be deduced in a straightforward manner from these algebraic 
identities. The local Hamiltonians of the XXZ spin chain satisfy the same algebra and the state space of each site 
is also two-dimensional: one goes from the modified matrix W to the XXZ spin chain with complex non-hermitian 
boundary fields through a simple change of variable discussed in details in [l^. For this reason, many methods 
introduced for example in [HI, [H, |29( for the XXZ spin chain can be applied successfully flGl , [TtI . to the ASEP. 
The next sections will focus however on the point of view and the language of stochastic exclusion processes. 



B. The stationary measure and the matrix ansatz 

The first eigenvector of the Markov transition matrix W that has been known \d\ is the stationary measure (eigen- 
value A = 0) when one discards any information about the current (s — 0). A configuration C of the system is a vector 
C = T1T2 . ■ - tl in {0, 1}^. The stationary probability of observing a configuration C is given by the matrix ansatz: 

1 ^ 

Pst{TlT2 ...Tl) = n (^^^ + (1 - ^')^) 1"^))' ("-11) 



i=l 



where the product of the matrices D and E is ordered from i = 1 on the left to i = L on the right and the vectors 
((yV| and |V)) are two vectors living in some auxiliary space, which is not the state space of the system. The condition 
VFjPst) = requires that the matrices D and E and the vectors {{W\ and |V)) satisfy the algebra Q: 

pDE-qED = D + E, (11.12) 
{{W\{-fD - aE) = ((>V|, (11.13) 
{f3D^6E)\V)) = |V)). (11.14) 

For generic values of the parameters, this algebra does not have any finite-dimensional representation, for which the 
vectors ((Wl and |V)) could have a physical interpretation. However, for some submanifolds of the parameter space 
@i @ , one can find finite-dimensional matrices D and E satis fying th e algebra above which have an interpretation in 
terms of superposition of shocks diffusing in the systems [13, ElllSl- I* would be interesting to understand the link 
with the eigenvectors obtained through the Bethe ansatz. 



C. Physical relevance of the eigenvectors for a non-zero parameter s 

The large deviation function f{j) of the current describes the probability of observing atypical values Q — jr of 
the current during a duration t. One way to understand these rare events is to characterize the time evolution of the 
system conditioned on the production of such an atypical current. One is interested in two particular characterizations 
of this conditioned evolution: the distribution of the final state knowing that an atypical current Q = jr was observed 
between and r on one side, and the conditioned transition matrix in the middle of the time interval [0, t] on the 
other side. 

The first case is easy and corresponds to the higher order terms of pi.9[) : the probability of observing a configuration 
C at T conditioned on j can be studied by a saddle-point analysis of the inverse Legendre transform of P(C, s\Co). 
One has 



- PriC,jT\Co) 



P^UC) = ^^T;.; ' :L. _ (C|m(^)) (11.15) 



where |/zi(s)) is the right eigenvector of W for s such that f'{j) — ~s. One sees thus that the physical interpretation 
of the first eigenvector of W is the description of the final state obtained after the observation of an atypical current 
j such that s = f'{j). 
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The second case can be obtained by cutting the time interval [0, r] into three intervals [0, ri], [ri,r2] and [T2,r]. 
One is thus interested in the probability -Pcond(C2, Q2, 721^1, Qi, ti; Cq, C, Q) that the system is in configuration C2 at 
T2 , with an integrated current Q2 , conditioned on the fact that the systems starts in Cq and ends in C with a current 
Q and knowing that the system was in Ci with a current Qi at the first intermediate time ri. This conditional 
probability can be written as the ratio of the probabilities of two histories, that can be further decomposed with 
Markov property: 



-Pcond(C2, Q2,T2\Ci,Qi,Ti;Co,C,Q) 



P(Co,0, O^Ci,Qi,ri ->C2,Q2,T2 -^C,Q,r) 
P(Co,0,O^Ci,Qi,ri ^C,Q,r) 

Pr-rAC,Q - Q2\C2)Pr,-rAC2,Q2 - Ql |Cl )P., (Ci , Qi [Cp) 
Pr-rdC,Q - Ql\Cl)PrdCl,Ql\Co) 

Pr-rAC, Q - Q2\C2)Pt,-tAC2, Q2 - Ql\Cl) 
Pr-rdC,Q-Ql\Cl) 



(11.16) 



where Pt{C, Q\C') is given in (|II.5[) . The dependence on Cq disappears because it is "erased" by the conditioning on 
having Ci at ti . If r2 — ti = Ar is fixed and ti and t — T2 go to infinity, the large deviation behaviour (jll.9[) shows 
that the dominant contribution to the conditioned probability is obtained when Qi and Q2 = Qi + AQ maximize 

(r - ri - At)/ - (t - n)/ 



Tl — At / \ T — Tl 

As expected Qi has to be of order jri, which means that the realization of the atypical current is distributed all over 
the time interval. A more careful study of the long time behaviour of (jll.9[) as for (jII.15|) shows that: 

Pcond(C2,Qi + AQ,ri + AT|Ci,Qi,ri;Co,C,jT) ~ M£lMe-^-Mi(.)+.AQp^^(^^^ ^qI^^) 

where s is given by f'{j) = — s and (/ii(s)| is the left eigenvector of the modified matrix W. The behaviour of the 

r.h.s. of the previous equation for an infinitesimal time Ar <C 1 shows that the transition matrix W^cond(^2,Ci) at an 
intermediate time conditioned on the production of a current j at a much larger time scale is given by : 

{fii[s)\Li) 

One checks that this new matrix has an eigenvalue A = that corresponds to the pseudo-stationary state conditioned 
to produce a current j. For j = jnwy, one has s ~ and thus /ii(s) = and {ni{s)\C2) does not depend on C2: one 
recovers the stationary measure of W given by the matrix ansatz. 

For s 7^ 0, the term (/j,i(s)|C2)/(/ii(s)|Ci) introduces an effective interaction between the particules that tends to 
enhance or reduce the jumping rates of particles, depending on their environment. A study of this interaction in a 
simple case is performed in section [V] 

More generally, the eigenvalues of W^^l^^ are given by Acond = A — ^i(s) where A is an eigenvalue of W for 
s = —f'{j)- The corresponding eigenvectors are given by: 

(C|Aeond) = (Mi(s)|C)(C|A) (11.19) 

where |A) is the right eigenvector of W for the eigenvalue A. 

This result shows that the characterization of the evolution of the ASEP conditioned to produce an atypical current 
j involves the right and left eigenvectors of the matrix W for a suitable value of the parameter s conjugated to the 
current. Moreover, the slowest relaxation times for this conditioned evolution is simply given by — (/X2(s) — ijli{s))~^ 
where /X2(s) is the second eigenvalue of W. One notices that the interpretation of the new matrix wj^^^ does not 
seem to have a simple physical counterpart in the quantum XXZ spin chain case. The next sections are devoted to 
the construction of some of the eigenvectors |A) by Bethe Ansatz methods. 



D. The coordinate Bethe Ansatz for periodic boundary conditions 



The second situation where the eigenvectors of the asymmetric exclusion process are known is the case of periodic 
boundary conditions. In this case, the number of particles N is constant and the 2^-dimensional state space splits 



6 



into L + 1 sectors of dimension L\/{N\{L — N)\) where N is the number of particles. Moreover, there is a duality 
between the sectors of N and L — N particles which corresponds to the particle-hole duality of the ASEP. Since there 
is no boundary, the current will be counted on each site, so that one must consider modified operators 



/O ON 

-q pe/ 

qe-' -p 

VO Oy 



(11.20) 



where s is the parameter conjugated to the current. In the framework of the Bethe ansatz, the state |00 . . . 0) is stable 
under all the operators Wi^i+i since 

m,^+i\o),\o),+i=o (11.21) 

and it can be used as a vacuum state on which particles (excitations) can be added. An isolated particle diffuses 
freely in the system with jumping rates p and q: 

w,,,+i|l)j|0),+i = pe'*|0),|l),+i -p|l),|0),+i, 
w,,,+i\0),\l),+i = (ze-^|l),|0),+i-g|0),|l),;+i. 

Thus, in the sector = 1, the eigenvectors are plane waves X^i^i ■ ■ ■ Oj-iliOi+i . . . Ol) with eigenvalues pe^ / z + 
qe~^z — {p + q), such that the parameter z satisfies z^ = 1 due to the periodicity of the lattice. The Bethe ansatz for 
general A^ consists (soj in assuming that eigenvectors are superpositions of plane waves with parameter z^ for each 
particle, such that an eigenvector \ip{zi, . . . , z^)) is given by: 

mz^,...,z^))= j2 E ^-<(i)<(2)---^:wio>r"'ii>.jo)^-i|i).,...io)^„+i. (11.22) 

X ord. {T(ESn 

The notations are defined as follows: x ~ {xi, X2, ■ ■ ■ , xn) is a vector ordered such that 1 < xi < X2 < ■ ■ ■ < x^ < L, 
the set Sn is the set of the A^! permutations of A^ elements and |0)^ stands for the tensor product |0)j ® . . . (g) |0)fc 
from site j to site k. 

By acting with W on configurations with no pair of neighboring particles, one checks that the eigenvalue 
A{zi, . . . , zjv) is given by the sum: 

A(zi,...,ZAr) =^ ( ^+9e-^zfe-p-g) . (11.23) 



fc=i 



When two particles are on adjacent sites, the exclusion effect breaks the simple diffusion of the particles and the 
coefficients have to be tuned so that \ip{zi, . . . , zn)) is still an eigenvector. For two particles with Bethe coefficients 
Zo-(fc) and ^^(fc+i) on sites i and i + I, one can allow the particles to exchange their Bethe numbers and, in order to 
have an eigenvector, one must have: 



+qe z„i^k)-p-q) A^z^f^^^z^^^^^^ /JJ24) 

Za{k+1) / ^ ' 



I p€^ „ \* 7 i I 1 

+ ( — + 96 z^(fe+i) -p-q) A„or^,,^,+^z^^^^-^)Z^^f^■^ 

^cr(fe) / 



ID^Ii) 



Since Wi^iJ^i\l)i\l)i+i = 0, one deduces immediately from the previous equation that the amplitudes must satisfy: 

^g-oTfc.fc + i _ Qpcriod (2^cr(fc+l), ^er(fc)) I'TT OK ^ 

2 — 7 \"! (^ii.^oaj 

^o- Operiod (^cr(fc) : 2<T(fc+l) } 

apcriod{z,z') = pe" + qe^'' zz' - {p + q)z. (II. 25b) 

If one starts from a permutation a and makes the first particle exchange its Bethe root ^^(i) with the second one, 
then makes the second particle exchange its Bethe root with the third one and so on until the last one, one obtains 
an amplitude ^ctoc2 ni where C2...N1 is a cycle of length L. Because of the periodic boundary conditions, the A'^-th 
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particle can also exchange its Bethe root with the first particle when they are on the sites L and L + 1 = 1 and one 
must recover the initial amplitude A„ up to a factor 2:^(1) ■ '^^^ Bethe roots Zj must then satisfy the equation for all 
j in {l,...,n}: 

z^^ = (-1)^-1 n (11.26) 

The resolution of these equations, at least numerically for finite N and L, or analytically for TV, L — > 00, gives 
the complete spectrum with (jII.23P : a detailed study of these equations (|II.26|1 was performed in [3l|, [s^l in relation 
with noisy Burgers equation and the six-vertex model. Eq. pi.22p gives the eigenvectors and allows one to compute 
correlation functions (see [H,!!!,!!!] for the XXZ spin chain). 

Equations similar to (jII.26P have been derived for special sets of parameters for the ASEP with open boundaries 
EE El- In this case however, the number of particles is not conserved at the boundaries and the integer N in 
(III.26| is replaced by some integer k that depends on the parameters a, /3, 7, (5, p, q and s. The procedure followed 
in [20I [21I [2^ is based on special prop erties of the transfer matrix and avoids the question of the eigenvectors. 
The algebraic Bethe ansatz used in (l9l . W\ for part of the spectrum of the XXZ spin chain was a first step to the 
determination of the eigenvectors. However the physical interpretation of the integer k in the context of the ASEP 
and simple approach such as the coordinate Bethe ansatz was still missing. In the next section, I show how one can 
determine with a simple procedure the excitations, which replace the particles, and how the integer k emerges in a 
straightforward way. In section llVi I tackle the problem of the scattering of these excitations at the boundary and 
the derivation of the Bethe equations. 



III. BULK INTEGRABILITY AND DESCRIPTION OF THE EXCITATIONS 



A. Product measure and the open chain 

For periodic boundary conditions, the particles |1) can be seen as excitations added on a reference vacuum state 
|00 . . . 0). In this case, the vacuum state is also an eigenvector (ground state) of W for the eigenvalue A = 0. For 
open boundaries at s = (the matrix W is stochastic), the ground state is given by the matrix ansatz (III.lip but, up 
to our knowledge, it has not been possible to use it as a vacuum state on which excitations may be added. Thus, we 
introduce a difference between the ground state which is the eigenvector whose eigenvalue has an extremal real part 
and the local vacuum states which are product states stable under the bulk operators and are used to separate 

excitations. 

If one looks for a tensor product |w)i|a;')i+i different from |0)i|0)i and |l)i|l)i+i such that 

Wi^i+i\uj)i\uj')i+i=Q, (III.l) 

then one sees that, up to multiplicative constants, the two vectors must be such that j39| : 

a |0), + C|1)„ 

« |0),+i+C(p/g)|l),+i, 

where C is still a free complex number. In the four-dimensional state space of two adjacent sites, it shows that a third 
tensor product satisfying pil.l[) exists besides |00) and |11). In the following sections, this number C will be tuned 
so that the vacuum states also behave well under the boundary operators Bi and B^. 
We thus introduce the vectors \uj{c))i on site i defined by : 

|c^(c)), = |0), + c(p/g)^|l),. (IIL2) 

These vectors satisfy by construction Wi,i+i\Lo{c))i\uj{c))i+i = and it will be useful to define the vacuum state from 
site i to site j through 

\n{c))l = \uj{c)), ® |w(c)),+i \lu{c))j. (III.3) 

One obtains consequently the stability of the vacuum state under the bulk dynamics: 

|C!(c))f = 0. (IIL4) 
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The boundary operators Bi and are still missing in (|III.4|) . If one wants the vacuum state [^(c))^ to be 
an eigenvector of W, then |r2(c))f must be an eigenvector of both Bi and B^. These two conditions determine 
independently twice the same coefficient c and give a constraint on the parameters. For each boundary operator, 
there are two choices for c, which correspond to the two possible eigenvectors. The sign of c(p/g)* is constant all along 
the chain and thus only two choices remain out of the four possibilities. The condition for the vacuum state |f2(c))f 
to be an eigenvector and the corresponding eigenvalue are thus summarized in the following table: 



Condition 


Eigenvalue 


Value of c 


(p/g)^-ie^ = 1 


A = 

A = -a-/3-7-(5 


c = e''{aq)/{-ip) 
c = -e^iq/p) 



(IIL5) 



One checks that these conditions are particular cases of the conditions obtained in [TtI. [TsI for which the number of 
Bethe roots in one of the two sectors is (A; = — L/2 or L/2 — 1 with the notations of U ll]). Thus the state pil.3|) 
corresponds to the case where there is no excitation and can play the role of a vacuum state as expected. 

In the case of the XXZ spin chain, the idea of finding vacuum states through local rotations already appeared in 
[sgI [37} and has been used most in the algebraic Bethe Ansatz framework [23, '37| . In the present case of the ASEP, it 
appears as a locally stationary two-sites state. In the case where c > 0, the state (|III.3|) can be seen as a densit y pr ofile 
with local Bernoulli measures of intensity c{p/q)^/{l + cijp/qY). In the language of shocks developed in llll IT^. 
it corresponds to L — 1 consecutive shocks. When c < 0, this interpretation in terms of density breaks down. 



B. Cutting the product measure and introducing excitations 

Excitations can be defined as local perturbations of the product measure (|III.3[) that diffuse freely under the bulk 
dynamics X]i=i^ when they are isolated. Thus, one can look for excitations \^l{c))i~^\4>)x\^{c))^^^i such that 

the action of X^^Li^ Wi,i+i on it gives a linear combination of |r2(c))^ x'\^{c))^i where x' — x ^ e where e = 
or ±1. However, no solution with this form can be found. 

A way of relaxing one of the constraint is to consider a state \VL{c))1^'^\(I))x\^{c'))x+i where c' ^ c. A solution is 
found if c' = c{q/p). In this case, the excitation \<j>)i takes the form: 

\^)^ = \0)^ + lyip/qy\^)^, (III.6) 

and its dynamics is the same as the one of a single particle among empty sites : it jumps on site i + 1 with rate 
p and on site i — 1 with rate q. In order to make the computations of section Hvl easier, it is useful to expand \(f))i on 
the two vectors \uj{c))i and \uj{c'))i (they form a basis 40j as soon a.s p ^ q) and replace the parameter which is still 
free by a parameter t, also independent of i, such that : 

m,c,c%=(^^^ {t\Lu(c)), + {l-t)\u;{c%). (111.7) 

where, once again, c' — {q/p)c. The factor (■\fqJpY is also introduced for later convenience such that the Bethe root 
of an excitation that has an initial Bethe root z and that is refiected becomes without additional factor even with 
p ^ q. The case t = makes the first vacuum state end at site i — 1 and the second start at site i] the choice t = 1 
makes the first vacuum end at site i and the next one start at site z + 1. Intermediate values of t give a superposition 
of both and seem redundant: we will see however in section [IV] that it can be useful to tune the value of t. 

The action of a single operator on an excitation involves the action of this operator on a tensor product 

\uj{c))i\u}{c'))i+i and one verifies that: 

■Wi^i+i\uj{c))i\uj{c'))i+i = p\uj{c'))i\uj[c'))i+i + q\uj{c))i\uj{c))i^i - [p + q)\uj{c))i\uj{c'))i+i. (III.8) 

One deduces that an excitation moves according to: 

Wi^i+i\uj{c))i\(j){t,c,c'))i+i = ./pq\cf){t,c,c'))i\uj{c'))i+i - q\uj{c))i\(j)[t,c,c'))i+i + \uj{c))i\v)i+i, (III. 9a) 
Wi,i+i\(j){t,c,c'))i\uj{c'))i+i = .Jpq\uj{c))i\(j){t,c,c'))i+i - p|0(t, c, c'))i|w(c'))i+i - \v)i\uj{c))i+i, (III.9b) 

where \v)i — {^J q/pYiq — p)\0)i. The opposite signs of the telescopic terms ±|u)i imply that they disappear up to 
boundary terms under the action of 
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The action of X^i^/ li'i.i+i 



on a state 



|f7(ci))r-i|0(ii, ci, C2))., |r!(c2)):j;}|0(i2, C2, C3))., . . . | r!(c„+i )) , 

with Xfc + 1 < Xfc+i (no excitations on adjacent sites) and Ck+i = {q/p)ck, thus makes each excitation at site Xi jump 
on the sites ± i or stay on the same site. Thus, the superposition of plane waves 



(III.IO) 



X ord. (tG»Sti \k—l / 

behaves, up to boundary terms^ as eigenvectors of the bulk dynamics 




= A, 



bulk 



t) + boundary terms, 



with the eigenvalue given by 



Abulk(2l, ■ 



■ T Zn) 



k=l 



Zk 



p- q 



= ^A(zfc), 

k=l 



(III.ll) 



(III.12) 



as long as the scattering of excitations on adjacent sites satisfies integrability conditions. The discussion of the 
boundary terms in (jlll.lip is the object of section ITVl and the discussion of the scattering of two adjacent excitations 
is performed in section UlI CI below. 



Several remarks have to be noticed at this stage. First, the vacuum states |^^(cfe))i^^_j^+i and \fl{ck+i)) 
the left or the right of an excitation are not the same since Ck+i = {q/p)ck =/= Ck, except for the symmetric simple 
exclusion process (SSEP). This difference explains why none of the eigenvectors can be used as a vacuum state and 
why we have introduced a distinction between the ground state and vacuum states at the beginning of section Hill 

The second remark relies on the identity |w(c„))i = |a;(c„+i))i+i. For t — or t = 1 (for other values of t, the 
excitations as in piL7p can always be decomposed), one sees that, near the site of the excitations, the states are lo cally 
products of two identical Bernoulli measures and one recognizes a similar structure as the one described in [l^, fTlf . 
The authors of [HI describe the states as combinations of shocks that separate Bernoulli product measure; here, 
on the contrary, the vacuum states are made of shocks separated by excitations that can be thought as Bernoulli 
product measures. This duality shocks/Bernoulli product measures seems to play a role similar to the particle-hole 
duality present in the periodic lattice when the number of particles is conserved. 

One can now attempt to construct an eigenvector from the state (jIII.10[) . If one discards the scattering of an 
excitation on one boundary (see section Hv]) and one requires that the first vacuum state \ fl{ci))i^~^ and the last one 
\rt{cn+i))x +1 eigenvectors of the boundary operators i?i and Bl, then the results (jlll.sp have to be replaced for 
n excitations by: 



Name 


Condition 


Eigenvalue A 


Value of ci 


Value of Cn+i 


(A) 


7(5 V?/ 


Abulk(^^l, • ■ • , Zn) 


ci = e''(ag)/(7p) 


cn+i = (5g^)/(/3p^) 


(B) 




Abuik(2:i, . . . , 2„) - (a + /? + 7 + <5) 


ci = -e^iq/p) 


Cn+1 = -{q^/p^) 



(III.13) 



The conditions presented in this table are exactly the ones obtained in [l^, [3, for which the spectrum is 
parameterized by Bethe roots Zk- The next sections are devoted to the determination of the Bethe equations, of the 
value of the U contained in the excitations and to the study of the boundary terms in (jllLlip . The final form of the 
coordinate Bethe ansatz is given in section HVl 



C. Scattering of two excitations on adjacent sites 



To have an eigenvector piLlip of the bulk dynamics up to boundary terms, one must check that the scattering of 
two excitations on adjacent sites \(j>{tk,Ck,Ck+i))i\(t>{tk+i,Ck+i,Ck+2))i+i is compatible with the isolated dynamics of 
the excitations. As for the case of periodic boundary conditions, the amplitudes A^- can be adjusted to satisfy this 
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log(ci(p/g)' 




los(H(H)^-i-"c=)=log(|) 



> site 



Figure 2: Schematic plot of the component (in logarithmic scale) along the basis vector of each vector in the tensor product 
\fl{ci))^-'-~^\(l){t,ci,C2))xi\^{c2))^^^\ ■ ■ ■ \^{cn+i))xn+i L vectors for n = 3: the component at site i is given by ci(p/g)*^™'' 
where rrii counts the number of excitations between sites 1 and i. The picture represents the case Ci > and t = with p > q. 



constraint, as in (jll.24p . However, from (jIII.9|) . one sees that additional terms \v)i and have to be introduced. 

Amplitudes and coefficients tk also have to be adjusted such that: 

W,.,+ i (^a4(fc)^^|fc+l) + ^<Torfc,, + i4(fc+l)^^tfc)) \4>{tk,Ck, Ck+l))t\(l>{tk+l,Ck+l,Ck+2))t+l 



1 



1 



,1+1 



icrorfc,fc + i^CT(fe+l)^<T(fc) 



(III. 14) 



A<yZ„(^k)^a(k+l) 



\4>{ik, Ck, Ck+l)) i\(l){tk+l, Ck+l, Ck+2)) i+l 
{\4>{tk,Ck,Ck+l))i\v)i+i - \v)t\(j){tk+l, Ck+1, Ck+2))i+lj ■ 



This system of four linear equations has non trivial solutions if and only if: 

tk — tk+i- 



(III.15) 



Thus, all the excitations along the lattice are characterized by the same global parameter t, which is still free. 
Moreover, one checks easily that the two amplitudes have to satisfy the same type of equation as (|II.25P for the 
periodic lattice, up to an irrelevant different normalization of the z^'s: 



A. 



a(z, z') 



O- (z<T(fc+l)7 ^CT(fc)) 

a (z<j(fc), ^o-(fc+i)) 

'pq + ^/Pqzz' - {p + q)z 



(III. 16a) 
(111.16b) 



These ratios relate the different amplitudes Aa- . To obtain the Bcthe equations in the periodic geometry, one starts 
with a given Aa-, permutes a Bethe root with all the other ones and then uses the periodic boundary conditions to 
recover the initial amplitude: the consistency conditions give the Bethe equations. In the present case, (IIIL16P allows 
one to permute Bethe roots and make one of them go from the first excitation to the last one. The scattering of an 
excitation on a boundary and the computation of its reflection coefficient are discussed in the next section. 



IV. SCATTERING OF THE EXCITATIONS AT THE BOUNDARIES 



A. Reservoirs and integrability 



The integrability of a quantum Hamiltonian or a stochastic transition matrix W can be seen as the knowledge of 
a non-trivial one-parameter family of matrices t{z) that commute with each other and contain the matrix W. The 
commutation relations [t(z), t(z')] = imply that they can all be diagonalized in a common basis and thus the matrix 
W itself also. However, there is no general procedure to construct these eigenvectors. In some cases as the ASEP on 
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the periodic geometry and the XXZ spin chain, the structure of the matrices t(z) provides creation and annihilation 
operators and a vacuum state from which one builds the eigenvectors with n particles or excitations from the ones 
with only n — 1 particles. 

For periodic systems, the systematic construction of the family t(z) relies on the algebraic properties of the local 
matrices Wi^i+i: one needs to find a family of matrices -Ri.i+i(z) that contains Wi^i+i and satisfies the so-called Yang- 
Baxter equations. These equations are conditions on the interaction between three bodies (three sites for the ASEP). 
For the XXZ spin chain or the ASEP, they can be deduced from the Temperley-Lieb algebra (III.10[) satisfied by the 

VJi,i+l- 

For the open chain, the systematic construction of the family t(z) was originally performed by Sklyanin [3^. The 
standard Yang-Baxter equations describe the integrability of the bulk dynamics. The boundary interaction operators 
Bi and Bj^ have to be integrable also: a family of matrices Ki{z) (resp. K]^{z)) is associated to each boundary 
site, contains the operator Bi and B^ and must satisfy reflection Yang-Baxter equations involving both the R and 
the K matrices 38]. In the case of the ASEP, the algebraic construction of Ki{z) and Kl{z) also relies on the 
Temperley-Lieb algebra (jll.lOp . 

However, the diagonalization in the open case of the family t(z) through a creation-annihilation algebra as for the 
ASEP with a periodic geometry is not possible in general. A creation-annihilation algebra and a vacuum state was 
found in the XXZ spin chain with non-diagonal boundary terms only at exceptional points in the parameter space . 
These exceptional points are exactly the same as in (jIII.13IIV.31l) . Attempts to avoid the creation-annihilation algebra 
that is not valid outside these points use other algebraic properties of the matrices t(u) (fusion rules, etc). However, 
they often lead to equations satisfied directly by the eigenvalues and prevent the construction of the eigenvectors 
[20I [21I [2^ . Other constructions for the XXZ spin chain have also been developed recently [2^ and may give also the 
full spectrum in terms of Bethe roots [1^ . 

The point of view followed here tries to avoid as far as possible the special algebraic properties of the model, 
so that the approach may be adapted more easily to other integrable models which do not necessarily rely on the 
Temperley-Lieb algebra. It may also help to identify new types of integrable boundary interactions for models whose 
bulk dynamics is already known to be integrable. 

Section [TV BI is devoted to the scattering of the first (resp. the last) excitation |(/)(t, ci, C2))i (resp. |(/)(t, c„, c„+i))i) 
on the left (resp. right) boundary; reflection and transmission coefficients are computed. Section IIV CI uses the 
expression of the reflection coefficients on the boundary to establish the Bethe equations satisfied by the roots Zk- 
The complete expression of the eigenvectors, i.e. the coordinate Bethe Ansatz, is presented in Section IIVDI and 
its consistency is also checked. Finally, the formalism is extended to left eigenvectors of W as well and the whole 
spectrum, as obtained in [l^ [l3, [3 , is described. 



B. Reflection of a first excitation on a boundary 

Sections nil Bl and nil CI have presented how excitations diffuse under the bulk dynamics. This section presents how 
an excitation gets scattered when it reaches a boundary. Computations are presented in details for the left reservoir 
but are valid for both boundaries (by changing a<->(5, 7<->/3, p<->q and x L -\- 1 — x). 

For closed boundaries {Bi = Bl = 0), the number of particles is conserved. In this case, coordinate Bethe ansatz 
takes a form similar to pi.22[) except that each Bethe root Zk can appear also with its inverse corresponding to 
a plane wave that propagates in the reverse direction. The amplitudes Aa- have to be replaced by amplitudes A^,? 
where e = (ei, . . . , ejv) is a vector of the hypercube with efe = ±1 representing the direction of propagation of each 
Bethe root. When a particle with Bethe root Zk arrives at site 1 or L, the boundary conditions induce a relation 
between Aa-^e and Aa-^r^e where rfce is the same vector as e* except that the fc-th component has a fiipped sign. The 
consistency conditions of all the amplitudes give the Bethe equations. 

In the case of the ASEP, the boundaries are not closed for the excitations when the parameters take generic values. 
Tuning the value of t introduced in l|III.7|) can close a boundary for the excitations (see below) but not the second 
one. Since there is no symmetry-related reason to close one or the other by tuning t, this parameter will be kept as a 
generic parameter. 

Since the boundary are not closed for excitations, an excitation that carries a Bethe root z and arrives at the 
boundary can either be reflected or removed. In the first case, it gives a second plane wave with Bethe root z~^ 
and an amplitude multiplied by a reflection coefficient Ri{z). In the second case, the state of site 1 is the second 
vacuum state |a;(c2))i with an amplitude multiplied by a transmission coefficient Ti{z). Moreover, in the case of the 
removal of an excitation, there are only n — 1 excitations that contribute to the bulk part of the eigenvalue piI.12[) : 
the determination of Ti{z) must take into account that Bi acting the second vacuum state |u;(c2))i must recover the 
missing contribution A(z) ^/pq{z + ^/z) — p — q to the bulk part of the eigenvalue. Moreover, depending on the 
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condition chosen in piI.13|) . the first vacuum state |ti;(ci))i is an eigenvector of Bi with eigenvalue Ai = or —a — 7: 
the action of Bi on the excitation or on the second vacuum state must also recover this contribution Ai . 

The boundary operator Bi induces a coupling between the three states Aiz\(j)(t,ci^C2))i, which corresponds to a 
plane wave with Bethe root z, AiZ~^\(/){t,ci,C2))i, which corresponds to a reflected plane wave with Bethe root 
and A[ = Ri{z)Ai, and the state A'-[\u} (02)) 1, which corresponds to the second vacuum state with A'( = Ti{z)A[. To 
have an eigenvector, the couplings must satisfy: 



Aiz\,p{t, ci, C2))i + A[z-^\<P{t, ci, C2))i + A'lXca))] 



= (Ai + ^-q^Aiz\(t){t,ci,c2))i + (^Ai + y^^z- q'^A'^z ci, C2))i 



(IV.l) 



+ (Aiz + A[z-^) \v)i + (Ai + A(z))A'iV(c2))j 



The first two terms of the r.h.s. are the contribution to the eigenvalue that is complementary to the one given by the 
action of u'1.2. The third term corresponds to the complementary part to the telescopic term \v)i left by wi^2 (see 
eq. pil.9p ). The fourth term is the contribution to the bulk eigenvalue that the second vacuum state must contain to 
compensate the disappearance of the first excitation. 

To solve pV.ll) . one must consider the action of Bi on the two vacuum states |ci;(ci))i and \uj{c2))i- The two cases 
presented in (|III.13p for which Ai = or —a — 7 can be treated simultaneously: 



Bi\uj{ci))i Ai|t^(ci))i 

B,Hc2)h = (-a - 7 - A,)Hc2)h + + ^^^^ + ^"'^ + k(ci))i 



(IV.2) 
(IV.3) 



Expanding (jlV.ip in the basis (|w(ci))i, |aj(c2))i) gives 



Ci 



(jp + q)hi + {aq + 7p) 



= P 



t + q{l-t) 



Az+p 



A(z) + a + 7 + 2Ai 



pqz{l — t) + pt 



Az- 



Mzt + q{l-t) A[z-^ 
(1 - t) +pt\ A[z~^ 



Ci 



The resolution of this system gives the value of A and A' as a function of the global normalization constant Ci 



Ci = (Aiz + A[z-^){l-t) + y^A'l 



pq z 



Hi 

q z 



+ (l-t) A^z 



A(z) + a + 7 + 2Ai 



[p + q)Ai + {aq + 7p) 



pqz 



7}j(^^tz + (l-t)^ A[z-^ = (^A(z)+a + 7 + 



2Ai 



(p + q)Ki + [aq + 7p) 



pq 



Up to a redefinition of the normalization Ci , one obtains : 

Aiz 
A\z-^ 
Vi{z) 



CiVi{z)pi{z) [z- 



CiVi{l/z)pi{l/z) 



A(z) + a + 7 + 2Ai 



1 



{p + g)Ai + {aq + jp) 



pqz 



J'-tz+{l-t) 

V 1 

The reflection coefficient of the plane wave with Bethe root z is given finally by: 

o . , ^ ^ _£Ml/£)pi(l/£) 
' A, z-W^{z)p,{z) 



(IV.4) 
(IV.5) 

(IV.6) 
(IV.7) 
(IV.8) 
(IV.9) 

(IV. 10) 
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Moreover, the computation of A'( as a function of Ci gives: 

2 



1 



p 



1 



pH^ + (1 - t)2 {pq -{p + 9)Ai - (ag + -fp)) + pt{l - t){p + q ^ a - - 2Ai) 

2 



ql-t 



(IV. 11) 



The second representation of the coefficient is ill-defined for t = or 1 since Vi{z) diverges as z 
A'l for t = or 1 is obtained by taking the limit i ^ or 1 in this second expression. 

On the second boundary, the three amplitudes have to satisfy an equation similar to (jIV.l|) : 



0. The value of 



Br 



A'l^z -^|0(i, c„,c„+i))l 



Alz |0(t, c„,c„+i))l 

Al + ^Jpqz - p)ALZ^\(j){t,Cn,Cn+l)) 



+ Al\L0{c„))l 



- {Alz^^ + A'^z-^) \v)i + {Al + Xiz))A'[\Loic„))L 
Similar computations as for the first boundary then give: 



Alz' 



A'^z-^ 



A'L 



Pl{z) 



ClVl{IIz)pl{z) ( z 
ClVl{z)pl{1/z) 

L . ^ 

z 

z 

fll 

J p z 



A(z)+/3 + 



t(l 

P 



t)]VL 



\l p 

{6p + I3q) 



t 



Cl 



pqz 



(IV. 13) 
(IV. 14) 

(IV. 15) 
(IV. 16) 
(IV. 17) 



From these expressions, computing the reflexion coefficient Rl{z) — A'j^/Al and the transmission coefhcient Tl{z) — 
A'I/Al is easy. The next section HV CI shows how to extract the Bethe equations from these reflection coefficient Ri{z) 
and Rl{z) and section llV DI gives the detailed form of the eigenvectors obtained from these reflexion and transmission 
coefficients. 

One can notice that the dependence on z of the transmitted amplitudes (jlV.lip and (jlV.lSP is all contained in the 
term {z — 1/^)^. The functions Vi and Vl depend only on the boundary rates a, (3, 7 and 5 and on the jumping 
rates p and q but are independent from the parameter t introduced in the definition of the excitations pil.7p . On the 
contrary, the functions pi and pL depend only on t and on the jumping rates p and q but are independent from the 
boundary rates. Thus the reflexion coefficients Ri{z) and Rl{z) have the form of a product of two terms, one that 
depends on the boundary rates but not on t and one that has the reverse dependence. 

One can also remark that the parameter s conjugated to the current is absent from all these expressions and appears 
only in the eigenvector through the coefficient that characterizes the vacuum states. 



C. Bethe equations 



If the transmitted terms left when an excitation leaves the system at one boundary are discarded in this subsection, 
the coordinate Bethe Ansatz for n excitations (when one of the condition piI.13p is satisfied) is given by 

n 

i^(.„...,z„))=_^ ^ _E<m4(ir 

X ord. crGtS„ e^Tin k—1 (^IV 18) 

X |f)(ci))r-'|0(i,Cl,C2)).J17(c2))^J;}|0(i,C2,C3)).J17(c3))^^;}...|f^(c„+i))^„+l 

+ terms with m < n excitations 
where 5„ is the group of permutations of n elements and Tin — {—1, 1}" is the n-dimensional hypercube. 
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If j = cr(l) and ej = +1, one can permute the Bethe root Zj of the first excitation with the second one ^^"2^)' 

and then with the third one and so on, until it reaches the n-th excitation: each time the two amphtudes v4^"'' 

and A*'"^ - are related through the ratio IIIII.I6I) . The n-th excitation, with Bethe root z, is reflected with a 

coefficient RL{zj) and the Bethe root becomes zj^ . One can permute this Bethe root in the reverse order with 
all the other ones until it reaches again the first one, and each permutation yields a scattering factor piI.16p . A 
reflection of the first excitation with Bethe root zj^ on the first boundary gives back the initial Zj up to a reflection 

coefficient Ri{z~^) = Ri{zj)~^. At this stage, one must recover the first amplitude ^o-.e and the cumulated product 
of all the scattering factors and the reflexion coefficients must be equal to 1. From the cycle of transformations 



> 

nil. 161 



A 

^CTOTi2,r„(i)e 



one obtains the Bethe equations: 



mil. 161 

> 

am. 161 



mil. 161 

> 

mil. 161 



A 



<TOC23.. 
i 



A 



O'0C23...7il,''<,(l)f 



(-1)"-! n 



a{zk,Zj) 
a{zj,Zk) 



n 



a{z, ^,Zk) 



(IV. 19) 



The contributions that comes from pi{z) and pl{z) cancel and the simplified Bethe equations are independent from 
the parameter t as expected: 



2L Vi{z,)Vl{z,) 
V,{z-')Vl(z-') 



n 



a{zj,Zk) a{zk,z^ ^) 
i{zk,Zj) a{z~^,Zk) 



(IV.20) 



This set of equations for j = 1, . . . , n is exactly the one obtained in [Tj, [3 ^^^y '^^^ P^^^ of the spectrum, after 
the change of variable z (zy^q/p— 1) / {^J q/p— z). If the first condition in (jIII.13|l is satisfied, one has Ai = A^ = 
and one recovers the second part of the spectrum in [l3] (eq. (2.13)). On the contrary, if the second condition in 
PII.13[) is fulfilled, one has Ai = —a — 7 and A^ = —(3 — 6 and one recovers the first part of the spectrum in [l3| 
(eq. (2.10)). I explain in section HV El how to obtain the second part of the spectrum in both cases by using the same 
approach on left eigenvectors of the transition matrix W . 

One can check that the equations pV.20|) have the expected symmetries. For example, changing a Zk to z^^ 
corresponds to exchanging the roles of a plane wave and its reflected partner and should not change the Bethe 
equations. One can check indeed that the set of equations PV.20|) is indeed invariant under the "gauge" transformation 
Zk — zl!° where the signs efe = ±1 are independent. 



D. Complete coordinate Bethe Ansatz and validity condition 

Eq. pV.lSp gives the Bethe equations by considering that the n excitations remain in the system up to reflexion 
coefficients at the boundaries. However section HVBI has shown that for generic values of t one must also consider the 
case where an excitation disappears at the boundary. These terms have to be added to the coordinate Bethe Ansatz 
to have correct eigenvectors. 

It is easy to see under which condition the boundaries are closed for the excitations: the two coefficients pV.ll|) 
and pV.lSP must vanish simultaneously. Introducing the reduced variable u = —^p/qt/{l — t) gives the condition: 

Vi{u)^Vl{1/u) = Q (IV.21) 

It is always possible to choose t and thus u so that the first term or the second term vanishes. However, they vanish 
simultaneously only under the additional assumption on the parameters that Vi{z) and Vi(l/z) have a common zero 
u* . For u = u*, there is no additional term in pV.lSp . 

For generic values of t, the additional terms in pV.lSp contains only n — 1 excitations that are moving in the bulk. 
In order to make computations easier, the free parameter t can be tuned to the value t = such that pV.lsp vanishes 
and excitations can leave the system only through the left reservoir. After the disappearance of \(j){t, ci, C2))i, the 
second excitation |(/)(t, C2, cs))^ can reach also the left boundary. To study its scattering, let us start with n excitations 
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as in PV.ISP and consider for simplicity the components a{j) = j for j > 3 and ej = +1 for j > 3, i.e. cr = id (identity) 
or cr = Ti2 (transposition of the the first two). There are 8 components corresponding to the possible permutations, 

ei = ±1 and 62 ~ ±1. Notations can be shortened by introducing the eight amplitudes A^g*^^ = A[2\ie2+ + 

^21 — ^Tl2,el£2+...+• 
Amplitudes >li2'^ and A2\^^ are coupled by piI.16[) . Moreover A^^^*^^ and ^12'''^^ are coupled by the reflexion 
coefficient i?i(zi) and both are coupled with the component, with amplitude ^[j*]27 i^i which \(j){t, ci, 02))^ has left the 

system with the root zi. In the same manner Aji^"'''' and ^21*' are coupled with Ri{z2) and are also coupled with 
the component, with amplitude ^p]!' corresponding to n — 1 excitations without Z2- Thus the four amplitudes A'^^lp 
and ^[2]i proportional to a global normalization constant K with: 

All^^z'^' = e2V^i(z^')pi(z^')a(zi,z5^)a(zf\zf)z2~^^ (^zi-l^if (IV.22) 

-eiVi{zl')pi{zl')a{z2,zl')a{z2\zl')z-'' (^2 - ^"j K (IV.23) 



^[2]l"^l 



If one removes then the second excitation, then the four amplitudes ^[j^]2 and ^p]! coupled to the same state, 
with amplitude ^[12] j which contains only n — 2 excitations with Bethe roots z^,. . . ,z„. The coupling of these five 
amplitudes is of the same type as in IjlV.ip except that one must recover the contribution of the Bethe root that has 
already disappeared from the system at the first stage: 



(IV. 24) 



Bi (^Aj^j^Zi + A[2]i2i ^ + A[*^]2Z2 + Aj^]2Z2 ^ j 1 0(i , C2 , C3)) 1 + A[i2] |w(c3)) 1 

(Ai + A(z2) + ^ - qjAf^^^zi + (Ai + A(z2) + VW^i - <i)a\2]iZi^ 

+ (Ai + A(zi) + ^ - 9)^[*1]2^2 + (Ai + A(zi) + VW22 - g) A[-^]222~^] \(l}{t, C2, C3))l 
+ (A+ ^Zi + Aj^jiZf 1 + 2^2 + ^fi]2^2"') I«>1 + (Ai + ^(^1) + A(Z2)) A[i2]|c^(c3))l 

The integrability of the boundary interactions with the reservoirs appears only at this point. Indeed, Eq. (jlV.ip 
can be written as three two-dimensional vectors, with three unknown amplitudes Ai, A'^ and A", whose sum is 0: 

Ai\a) + A[\a') + A'l\a") = (IV.25) 

In the generic case, |a') and \a") are independent and it is always possible to find A[ and A'( to satisfy this condition. 
On the contrary, pV.24[) can be written as a sum of five two-dimensional vectors that should give but four amplitudes 
are already fixed (up to an overall normalization constant) by the previous step. The sum can be written formally: 

+ A[2]ll^2"l) +4]2l&i2> +^ri]2l^r2>) + Al2]|fc) = (IV.26) 

where the vector formed by the first four vectors is fixed by the previous step. A solution for ^[12] exists if and only 
if \b) is proportional to the sum of the four other vectors. If it is not the case, the Ansatz breaks down and the Bethe 
equations pV.20[) are meaningless. For the boundary operator Bi and Bl defined in pi.2p . a lengthy computation 
shows that the sum of the first four vector is proportional to |6) and a solution exists indeed for A[x2]- 

For boundary operators, the integrability condition on the scattering of two excitations in the bulk piI.14p is 
replaced by a condition on the transmission coefficients of the first excitations zi and Z2' eq. pV.24[) states that they 
must be compatible with the reflection coefficients of the second excitation with the same two Bethe roots. The same 
procedure holds for the removal of the third excitation and, from integrability, one expects it be true until the removal 
of the last excitation \(t>(t, c„, c„+i))i although showing it rigorously would require more efforts. 

Thus for t = t*j^ such that the transmission coefficient pV.15l) at the right boundary is 0, the complete coordinate 
Bethe ansatz can be written as: 

n n 

i^(zr,...,z„))=5: E E E n 4(ir 

m = ^(n-m) aeSn-m (GHn-m k=m + l 

X |f](c„i+i))i" + '"V(^7Cm-|-l,Cm+2))2;™+i|f^(c,n+2))!^"+i"^|(/'(t, Cm+2,C„i+3))a;„,+2 . . . \Vl{Cn+l))x^+l 

(IV.27) 
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where m counts the number of excitations that have been removed at the left boundary. The (n — m)-dimensional 
vector ^("-™) — {xm+i, ■ ■ ■ ,Xn) with Xk < x^+i give the position of the n — m excitations that remain in the bulk. 
As before, the sets Sn-m and Tin-m are the group of permutations of n — m objects and the hypercube { — 1, 1}"^™. 

If t was chosen to be equal to i* so that the left boundary is closed instead of the right one, the coordinate Bethe 
ansatz (jIV.27|) would have been a combination of tensor product vectors starting with the vacuum |r2(ci))i and ending 
with the vacuum |f2(c„-m+i)) l- Both expressions should be equal up to a normalization constant. 



E. Second part of the spectrum and left eigenvectors 

In section IIV C[ only one part of the spectrum described in [l3, [13 is obtained with the ansatz (jlV.lSP . Up to 
now, it has not been possible to write down the ansatz that corresponds to the right eigenvectors in the second part 
of the spectrum. However it is possible to obtain the Bethe equations for this second part by considering the left 
eigenvectors of the matrix W. To simplify the notations, we will write them as right eigenvectors of the transposed 
matrix and the change |0) i— > (0| and |1) (1| allows one to go easily from one to the other. Up to a change 
of basis, W"^ with a given value of s is the same as W with a value s' = —rj — s where = [a(3) / {'-^5){p/ q)^~^ . 
This symmetry is known as the Gallavotti-Cohen symmetry and W'^ is related to the time-reversed properties of the 
exclusion process. 

From the definition 

= |0), + c|l)„ (IV.28) 

one checks easily that 

<.+i(p(5))^®P(5)).+i) =0 (IV.29) 
and the vectors \u){c))i can be used to build a vacuum state. Moreover excitations defined by 

|0(f,c,c')). - (^yi) {t\CJ{^)^ + (1 - imc'))^, (IV.30) 

where c' = {q/p)c, move with the same dynamics as in pil.9p : the jumping rates are the same and the telescopic term 
\v)i has the same expression. A state with n' excitations inside the bulk that separate vacuum states with coefficients 
Cfc characterized by Ck+i = {q/p)ck (as in pil.lOp ) is an eigenvector of the boundary operators I3f and Bj^ if one of 
two constraints below are satisfied: 



Name 


Condition 


Eigenvalue A 


Value of ci 


Value of Cn+i 


(A) 




Abuik(^i, . . . , Zn>) - (a + /3 + 7 + (5) 


ci — —6^*7/0; 


Cn+i = -(i/5 


(B) 




Abulk(2l, ■ ■ ■ ,Zn') 


Ci = e^^ 


2,1+1 = 1 



(IV.31) 



These constraints are exactly the same as in (|IIL13p except that L — 1 — n is replaced by n' . Thus, if the first constraint 
of (jIII.13[) is satisfied for given n then it is possible to build: 

1. right eigenvectors with n excitations of type (|III.7[) separating vacuum states \uj{ck))i and with eigenvalue 
Abuik(2:i, • ■ • , Zn) where the ^^'s satisfy the Bethe equations (|IV.20[) . 

2. left eigenvectors with n' = L — 1 — n excitations of type (|IV.30P separating vacuum states \Lb{ck))i and with 
eigenvalue Abuik(-Si, • ■ • , ^n') — (a + /3 + 7 + (5) where the z^'s satisfy a second set of Bethe equations similar to 
the ones that correspond to the other part of the spectrum in [l3, [HI . 

On the contrary, if the second constraint of piI.13[) is satisfied, the correspondence to the two parts of the spectrum 
in [ill is reversed. 

It is interesting to compare these results with the Matrix Ansatz pi. Ill) for s — 0. The construction of left and 
right eigenvectors described in this paper gives the complete spectrum described in jl^, [13, [3 t>ut only a subset of 
the right or the left eigenvectors. For s = 0, it is always possible to choose n = L — 1 within the constraint (B). The 
right eigenvectors describe the eigenvalues A(zi, . . . , zl-i) — (a + /? + 7 + (5), which do not contain the ground state 



17 



Generic case s 7^ 


Special case s — 


Part of the spectrum 


Left eigenvect. 


Right eigenvect. 


Part of the spectrum 


Left eigenvect. 


Right eigenvect. 


I (resp. II) 


CBA 


? 


A = 


(8)ti((o| + {i|) 


Matrix Ansatz 


II (resp. I) 


? 


CBA 


Other A's 


7 


CBA 



Table I: Summary of the eigenvectors described in this approach and the non-overlap with the Matrix Ansatz for s = 
(CBA=coordinate Bethe Ansatz described in this paper). In the first table, the part of the spectrum described by a set of 
eigenvectors depends on the constraint (|III.13P that is satisfied. The second table is obtained by taking n — L — 1 with constraint 
(B). 



A = described by the matrix ansatz pLlip . On the contrary, the construction explained in this paper give only one 
left eigenvector, which is a product state made of (0| + (1|, with eigenvalue A = 0. It corresponds to the conservation 
of the total probability for s = and is thus the ground state. A short summary of these results is presented in 
table jH. There is one known case [8, .9] for which the matrix ansatz does not give straighforwardly the stationary 
state and it corresponds precisely to the condition (A) in (jIII.13[) for s = 0: in this case, the coordinate Bethe ansatz 
corresponding to the constraint (A) in piI.lSp gives it with all the z^s going to 1. The eigenvectors presented here 
are never redundant with the Matrix Ansatz. 



V. APPLICATION: DYNAMICS CONDITIONED ON THE CURRENT IN THE CASE n' = 1 

The previous section shows how to build some of the left eigenvectors needed for the study of the dynamics W^^]^^ 
conditioned to produce a long time current j. Eq. pi.lSp shows that the first left eigenvector of W is needed. If this 
eigenvector belongs to the part I of the spectrum described in fig. HI then it is given by a coordinate Bethe Ansatz. 
For simplicity, only the case s = — hi{p/q), which corresponds to n = L — 2 and n' = 1, will be presented here: it is 
the first non-trivial case after s = since some of the left eigenvectors depend only on one Bethe root, noted z from 
now on. The corresponding left eigenvectors are obtained from pV.28IIV.30p and are given by: 



z^ + A'z~n I (K)(0| + 5(l| I ( m + m) (g) (0| + 5^(11 




(V.l) 



where c — p/q and i> is chosen in order to close the right boundary as in section HVDI The coefficients A, A' and A" 
are the ones given by the reflection condition at the left boundary in section HVBI The Bethe root z is a root of the 
polynomial equation 

^^iz)V^ . 1. (V.2) 
V^{1/z)Vl{1/z) ^ ' 

It is not trivial to see that the first eigenvalue /ii(s) of W belongs to this part of the spectrum, even for n' = 1. 
However numerical checks up to L = 6 show that it is the case: it will be assumed that it is the case for general L 
and z is the corresponding solution of the previous equation. For a configuration C = ti . . . where = or 1 , the 
component of the first eigenvectors of W on \C) is given by: 



(Mi(.)ic)=5:(A.^+A'.-) u/f ^^'( n i'l) n (v.3) 




This simple expression where all the numbers are known can be used directly to simulate the conditioned dynamics 
for the non-trivial value of s = — ln(p/g). 

To illustrate the effect of the conditioning, one can look at the modified jumping rate of a single particle in the 
system. For a configuration Cy = . . . 010 ... where the single particle is at position y, the particle can hop on the 
left side with probability qe^^U{y — 1)/U{y) or to the right with probability pe^^U{y + 1)/U{y). The contribution 
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U{y) = {fii{s)\Cy) acts as an exterior potential and is given by: 



Uiy) = {Ms)\Cy) = A" ( ) +J2{Az^ + A' z-^) 



cS 



y<x 



-Ik 



(V.4) 



Similar computations with two particles in the system gives the interaction, which may not be short-range, between 
two particles in addition to this exterior potential. More detailed results about this interaction will be presented in a 
future work. 



VI. CONCLUSION 



This paper shows how to construct eigenvectors of the transition matrix of the asymmetric exclusion process with 
different reservoirs at both boundaries modified to count the current, when the parameters satisfy some algebraic 
conditions. The Ansatz (|IV.27|) shows that one must allow excitations to leave the system through the boundary 
reservoirs: the integrability condition takes the form of a compatibility condition pV.26|) between the transmission 
coefficients of the m-th excitation and the reflexion coefhcients of the (m + l)-th one. 

A key feature in the identification of the special points specified by the constraints (jIII.13[) is the existence of a 
free parameter c in the local stationary measures \Lu{c))i that can be adjusted to the boundary operators Bi and Bj^. 
When c is positive, it is related to the mean density of particles on the sites of the lattice. An excitation can be seen 
as a moving frontier between two local stationary measures with different parameters c. The detection of the special 
points where a coordinate Bethe Ansatz (jIV.27|l is valid is then made easy once the relation between the left and 
parameters c and c' is known: one starts with ci fixed by the left boundary, computes C2,. . . ,c„+i with the previous 
relation and tries to match it with the right reservoir. 

The second step consists in studying how excitations move in the bulk and interact, with the usual technique of 
the coordinate Bethe Ansatz. At the boundary, the determination of the reflection and transmission coefficients for 
the first excitation is easy and it is enough to get the Bethe equations if one already knows that the coordinate Bethe 
Ansatz is valid. Checking its validity, and thus the integrability, requires however to check that the transmission 
coefficient of the first excitation are compatible with the reflection coefficient of the second one, etc. Once the Bethe 
equations are obtained, the ground state and the gap can be computed after some work [li,[l3,[l|. 

The advantage of the method presented here is that it does not use special algebraic properties of the model and 
one may hope that it could be extended to other integrable models. Moreover, the approach followed here provides 
an easy way of finding the submanifold of the parameter space where this coordinate Bethe Ansatz is valid; however 
this restriction to a submanifold appears only on presence of two boundaries: most of the construction may be valid 
without restriction for a half-line with only one reservoir. 

When the measure of the current is discarded by setting s = 0, interesting relations with the Matrix ansatz appear. 
The only right eigenvector that is not given by the coordinate Bethe Ansatz is the one given by the matrix ansatz. 
When constraints (A) in pil.isp is satisfied, then the matrix ansatz is known to fail in this case ([1, [^). Moreover, 
the decomposition of the eigenvectors into shocks and Bernoulli product measures is the one that appears when the 
matrix ansatz has finite representations. The parameter space contains submanifolds with a flnite dimensional matrix 
ansatz (finite number of shocks) and submanifolds (jIII.13[) with a finite number of Bethe excitations. From these 
considerations, the missing eigenvectors in the table |T] may be related to a modified form of the matrix ansatz and it 
would be interesting to know how the two approaches are related. Coupled with (|II.19[) , the knowledge of the missing 
eigenvectors would allow one to understand the full dynamics of ASEP conditioned to give a given current and may 
help to study the optimal density profile to produce this current. 

One may be interested in extending this construction to other models. The Bethe roots Zk are attached to the 
excitations in the eigenvectors and the spectrum is completely determined by them. The structure of the vacuum 
states does not play any role once the dynamics of their frontiers is known. It may be interesting to try to use other 
types of locally stationary measures as vacuum states in other models and check if their boundaries have integrable 
motions. 

A challenge is the determination of the structure of the eigenvectors for generic value of the parameters. The works 
[2^[25j are first steps in this direction for the XXZ spin chain. It would be interesting to find a physical interpretation 
in terms of shocks and Bernoulli product measures of these results in the case of the ASEP. It would also be very useful 
to see if the parametrization of the spectrum done [1^ in terms of Bethe roots also has a structure that generalizes 
the one presented here. 

Finally, it would be interesting to extend this formalism to other fields of statistical physics. Two active domains 
where integrability plays an important role and for which boundary conditions are relevant are tilings and loop 
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models. The knowledge of the eigenvectors and of the nature of the excitations in this models may provide interesting 
information about the typical configurations of these other models. 
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